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IN THE CLAIMS ; 

Please find a listing of the claims below. The statuses of the claims are shown in 
parentheses. 

1. (Canceled). 

2. (Previously presented) A method for gamut mapping of an input image using a 
space varying algorithm, comprising: 

receiving the input image; 

converting the color representations of an image pixel set to produce a corresponding 
electrical values set; 

applying the space varying algorithm to the electrical values set to produce a color- 
mapped value set; 

reconverting the color-mapped value set to an output image; and 

wherein the space varying algorithm minimizes the following variational problem: 

E{u) = \{D 2 + apD\ 2 )dQ, subject to we & wherein Q is a support of the input 
n 

image, 0 is the target gamut, a is a non-negative real number, D=g*(u-uo), g is a normalized 
Gaussian kernel with zero mean and a small variance a, Uo is the input image, and u is the 
output image. 
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3. (Original) The method of claim 2, further comprising: 
solving the variational problem at a high value of a; 
solving the variational problem at a low value of a; and 
averaging the solutions. 

4. (Original) The method of claim 3, wherein the step of averaging the solutions 
comprises using a spatially adaptive weighting scheme, comprising: 

Uf m ai[kJ] = w[A: 5 /*]w 5ma //[^ I /'](l-w[A: l /*])w/ 1/ ^[A: 5 /] ? 
wherein the weight w[kj] 9 comprises: 

w[kj] = — 1 -5- , and 

1 + /3\Vg*u 0 \ 

wherein P is a non-negative real number. 

5. (Original) The method of claim 2, wherein the variational problem is solved 
according to: 

du 



dt 



= ag * AD - g * D, subject to u e A 



6. (Original) The method claim 2, wherein the space varying algorithm is solved 
according to: 

u i/ l= u v + T ( aL V " sub J ect to u y G & > wherein 
dt, 

D n =g*g*(u"-uJ 
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L n = D 2 * (u"-u 0 ) and 

D 2 =g x *gx+gy *gy 

7. (Previously presented) A method for gamut mapping of an input image using a 
space varying algorithm, comprising: 

receiving the input image; 

converting the color representations of an image pixel set to produce a corresponding 
electrical values set; 

applying the space varying algorithm to the electrical values set to produce a color- 
mapped value set; 

reconverting the color-mapped value set to an output image; and 

wherein the space varying algorithm minimizes the following variational problem: 

E(u) = J (/j (D) + ccp 2 ( | VD | ))dO, subject to u e <9, wherein p x and p 2 are scalar 

functions, Q is a support of the image, & is the target gamut, a is a non-negative real number, 
D=g*(u-uo) ? g is a normalized Gaussian kernel with zero mean and a small variance a, uo is 
the input image, and u is the output image. 

8. (Original) The method of claim 2, further comprising: 

decimating the input image to create one or more resolution layers, wherein the one or 
more resolution layers comprises an image pyramid; and 

solving the variational problem for each of the one or more resolution layers. 
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9. (Previously presented) The method of claim 2, wherein the method is executed in 
at least one of a camera and a printer. 

10. (Previously presented) The method of claim 7, wherein the method is executed in 
at least one of a camera and a printer. 

Claims 11 and 12. (Canceled). 

13. (Previously presented) A computer-readable memory for color gamut mapping, 
comprising an instruction set for executing color gamut mapping steps, the steps, comprising: 

converting first colorimetric values of an original image to second colorimetric 
values, wherein output values are constrained within a gamut of the output device; using a 
space varying algorithm that solves an image difference problem; and 

optimizing a solution to the image difference problem, wherein the image difference 
problem is represented by: 

E(u) = \(D 2 + a\ VD\ 2 )dn 

Q 

subject to u g 9 , wherein Q is a support of an input image, a is a non-negative real number, 
S is the target gamut, D=g*(u-uo), g is a normalized Gaussian kernel with zero mean and 
small variance c, uo is the input image, and u is an output image. 
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14. (Previously presented) The computer-readable memory of claim 13, wherein the 
instruction set further comprises steps for: 

solving the image difference problem at a high value of a; 
solving the image difference problem at a low value of a; and 
averaging the solutions. 



15. (Original) The computer-readable memory of claim 14, wherein averaging the 
solutions comprises using a spatially adaptive weighting scheme, comprising: 

wherein the weight w[kj] 9 comprises: 

w[kj] = —± -5- , and 

1 + y#|Vg* u 0 \ 

wherein P is a non-negative real number. 



16. (Previously presented) A computer-readable memory for color gamut mapping, 
comprising an instruction set for executing color gamut mapping steps, the steps, comprising: 

converting first colorimetric values of an original image to second colorimetric 
values, wherein output values are constrained within a gamut of the output device; using a 
space varying algorithm that solves an image difference problem; and 

optimizing a solution to the image difference problem, wherein the image difference 
problem is represented by: 

E(u) = J (/J (D) + ap 2 ( | VD | ))<i/2, wherein p x and p 2 are scalar functions. 
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17. (Currently amended) A computer-readable memory for color gamut mapping, 
comprising an instruction set for executing color gamut mapping steps, the steps, comprising: 

converting first colorimetric values of an original image to second colorimetric 
values, wherein output values are constrained within a gamut of the output device; using a 
space varying algorithm that solves an imag e diff e r e nc e probl e m the following variational 
problem: 

E(u) = J(£> 2 + «|vZ)| 2 )rfQ, subject to we ,9, wherein Q is a support of the input 
n 

image, & is the target gamut, a is a non-negative real number. D=g*(u-un\ g is a normalized 
Gaussian kernel with zero mean and a small variance g, uq is the input image, and u is the 
output image ; and 

optimizing a solution to th e imag e diff e r e nc e probl e m, wherein the instruction set 
further comprises steps for: 

decimating the input image to create one or more resolution layers, wherein 
the one or more resolution layers comprise an image pyramid; and 

solving the imag e diff e r e nc e variational p roblem for each of the one or more 
resolution layers. 
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18. (Currently amended) The computer-readable memory of claim 1 7, wherein the 

instruction set further comprises steps for: 

(a) initializing a first resolution layer; 

(b) calculating a gradient G for the resolution layer, the gradient G comprising: 

G = A(w - Uq) + olk(u - uo\ wherein Ax is a convolution of each color plane of x with 



K LAP ~ 



0 1 0 

1 -2 1 
0 1 0 



and a k - a 0 * 2 2(k 0 , wherein is a Laplacian kernel, k is a 



kernel, and a 0 is an initial non-negative real number; 

(c) calculating an normaliz e d st ee p e st d e sc e nt iteration layer value Lj = Lj.i - |xo * 
Unsd * G, wherein |io is a constant , i is a specific resolution layer, and [iN gn is a normalized 
steepest descent parameter ; 

(d) projecting the value onto constraints Proj^Lj), wherein Proj^x) is a projection 
of x into a gamut and 

(e) for a subsequent resolution layer, upscaling from one resolution layer to another 
and repeating steps (b) - (d). 



19. (Currently amended) A method for image enhancement using gamut mapping, 
comprising: 

receiving a input image; 

from the input image, constructing an image pyramid having a plurality of resolution 

layers; 

processing each resolution layer, wherein the processing includes completing a 
gradient iteration, by: 
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calculating a gradient G , for G = A(u-un) + ak(u-un)« wherein u is the iterated 
image, up is the input image at the appropriate resolution layer, and a is a non-negative real 
number ; 

completing a gradient descent iteration; and 

projecting the completed gradient descent iteration onto constraints; and 
computing an output image using the processed resolution layers. 

20. (Canceled). 

21 . (Currently amended) The method of claim 19 5 wherein completing the gradient 
descent iteration (LJ) comprises calculating: 



wherein junsd is a normalized steepest descent parameter, is a constant, k is a number of 
resolution layers in the image pyramid, and j is a specific resolution layer. 

22. (Original) The method of claim 19, wherein projecting the completed gradient 
descent iteration onto the constraints is given by: 
Lj=Proj a (Lj), 

wherein Proj a (x) is a projection of x into a gamut S. 




Lj = Lj.\ - jj q • Mnsd' g > 
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23. (Currently amended) The method of claim 19, wherein constructing the image 
pyramid, comprises: 

smoothing the input image with a Gaussian kernel; 
decimating the input image; and 

setting initial conductiv e condition Lo = max {Sp} 5 wherein Sp is an image with the 
coarsest resolution layer for the image pyramid. 

24. (Original) The method of claim 23, wherein the Gaussian kernel, comprises: 

16 8 16 
Hi 
8 4 8 
HI 
16 8 16_ 

25. (Previously presented) The method of claim 19, wherein processing each 
resolution layer further comprises applying a space varying algorithm to minimize the 
following variational problem: 

4 

E{u) = J (D 2 + a | VD | 2 )<£f2, subject to w e 5, wherein Q is a support of the 

image, 0 is the target gamut, and D=g*(«-W0), wherein g is a normalized Gaussian kernel with 
zero mean and small variance a, uq is the input image, u is the output image, and wherein a is 
a non-negative real number. 



XT — 
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26. (Original) The method of claim 19, wherein processing each resolution layer 
comprises applying a space varying algorithm to minimize a variational problem represented 
by: 

E(u) = J (p (D) + apA I V£> | ))df2, subject to u e 9, wherein p and p 
are scalar function. 

27. (Original) The method of claim 26, wherein pi and p2 are chosen from the group 
comprising p (x) = \x\ and p (x) = Vl + x 2 . 
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